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Using the thermal field dynamics theory to convert the thermal state to a "pure" state in doubled 

f t 

Fock space, it is found that the average value of e-*" " under squeezed thermal state (STS) is 
just the generating function of Legendre polynomials, a remarkable result. Based on this point, 
the normalization and photon-number distributions of m-photon added (or subtracted) STS are 
conviently obtained as the Legendre polynomials. This new concise method can be expanded to the 
entangled case. 



I. INTRODUCTION 

Nonclassicality of optical fields is helpful in under- 
standing fundamentals of quantum optics and have many 
applications in quantum information processing [ij. To 
generate and manipulate various nonclassical optical 
fields, subtracting or adding photons from/to traditional 
quantum states or Gaussian states is proposed 043- 
For example, the photon addition and subtraction have 
been successfully demonstrated experimentally for prob- 
ing quantum commutation rules by Parigi et al. [7] . Re- 
cently, photon-added (-subtracted) Gaussian states have 
received more attention from both experimentalists and 
theoreticians since these states exhibit an abun- 

dant of nonclassical properties and may give access to 
a complete engineering of quantum states and to funda- 
mental quantum phenomena. 

Theoretically, the normalization factors of such quan- 
tum states are essential for studying their nonclassical 
properties. Very recent, Fan and Jiang [20| present a 
new concise approach for normalizing m-photon-added (- 
subtracted) squeezed vacuum state (pure state) by con- 
structing generating function. Hovifever, most systems 
are not isolated, but are immersed in a thermal reser- 
voir, thus it is often the case that we have no enough 
information to specify completely the state of a system. 
In such a situation, the system only can be described by 
mixed states, such as thermal states. In addition, the 
squeezed thermal states (STSs) can be considered as the 
generalized Gaussian states. 

In this paper, we shall extend this case to the mixed 
state, i.e., by using the thermal field dynamics (TFD) 
theory to convert the thermal state to a "pure" state in 
doubled Fock space, we present a new concise method for 
normalizing photon-added (-subtracted) squeezed ther- 
mal states (PASTSs, PSSTSs) and for deriving their 
photon-number distributions (PNDs) which have been 
a major topic of studies on quantum optics and quan- 
tum statistics. It is found that the normalization factors 
and PNDs are related to the Legendre polynomials in a 
compact form. 



Our paper is arranged as follows. In section 2, based on 
Takahashi-Umezawa TFD, we convert the thermal state 
Pth to a pure state in doubled Fock space by the par- 
tial trace, pth = tf [|0(/3)) (0(/3)|] (see Eq.® below). In 
section 3, we introduce the STS Si |0(/3)) with 5*1 being 
the single-mode squeezing operator for the real mode. 
It is shown that the average value of e^" " under STS 
is just the generating function of Legendre polynomials, 
a remarkable result. Based on this point, in sections 4 
and 5, the normalization factors and PNDs of m-photon 
added (or subtracted) STS are obtained as the Legendre 
polynomials, respectively. The last section is devoted to 
drawing a conclusion. 



II. REPRESENTATION OF THERMAL STATE 
IN DOUBLED FOCK SPACE 

We begin with briefly reviewing the properties of ther- 
mal state. For a single field mode with frequency a; in a 
thermal equilibrium state corresponding to absolute tem- 
perature T, the density operator is 



Pth = 
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where ric = [exp(fia;/(fcT)) — being the average 

photon number of the thermal state pth and k be- 
ing Bltzmann's constant. Note \n) — a'"'/Vri\\0) 
and the normally ordering form of vacuum projector 
|0) (0| — : exp{—a^a): (the symbol : : denotes the nor- 
mal ordering), one can put Eq.(IT]) into the following form 
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the last step, the operator 
exp [(e'^ — l) a^a] : is used. 



(2) 



identity 



where in 
exp (Aa^a) 

Recalling the thermal field dynamics (TFD) intro- 
duced by Takahashi and Umezawa (2ll - [23j |. its elemental 
spirit is to convert the calculations of ensemble averages 
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for a mixed state p, (A) = tr {Ap) /tr (p) , to equivalent 
expectation values with a pure state |0(/3)), i.e., 

(A) = (0(/3)|A|0(/3)), (3) 

where /3 = 1/kT, k is the Boltzmann constant. Thus, 
for the density operator pth, by using the partial trace 
method m, i.e., pth = [|0(/3)) (0(;9)|] where t? de- 
notes the trace operation over the environment freedom 
(denoted as operator a^), one can obtain the explicit ex- 
pression of |0(/3)) in doubled Fock space, 

|0(;3)) = seche'exp [a^at tanh^] |00) = S (6) |00) , (4) 

where |00) is annihilated by a and a, [a, a^] = 1, and S (6) 
is the thermal operator, S (9) = exp [6* (a^a^ — aa)] with 
a similar form to the a two-mode squeezing operator ex- 
cept for the tilde mode, and is a parameter related 
to the temperature by tanh^ = exp(— ^^). |0(/3)) is 
named thermal vacuum state. 

Let Tr denote the trace operation over both the system 
freedom (expressed by tr) and the environment freedom 
by tr, i.e., Tr = trtr, then we have 

tr{Ap„,)=Tr[A\Om (0(/3)|] 

= tr[A^\0m{0m], (5) 

and the average photon number of the thermal state pth 
is 



n, = Tr [a^alOm (0(/3)|] 



sinh^ 9. 



(6) 



Here we should emphasize that tr |0(/3)) (0(/3)| ^ 
(0(/3)| 0(/3)) , since |0(/3)) involves both real mode a and 
fictitious mode a. From Eqs.Q and ^ one can see that 
the worthwhile convenience in Eq. Q is at the expense of 
introducing a fictitious field (or called a tilde-conjugate 
field) in the extended Hilbert space, i.e., the original op- 
tical field state \n) in the Hilbert space H is accompanied 
by a tilde state \h) in H. A similar rule holds for oper- 
ators: every Bose annihilation operator a acting on H 
has an image d acting on H. These operators in H are 
commutative with those in H. 



where we have used 5'i (A) |0) =sech^/^Aexp[— a'l'^/2tanhA] |0) 
Further, note e^'^°'^ e~^°'°'^ — a^-|-ra^, and for operators 
A,B satisfying the conditions [A, [A,B]] = [B, [A, B]] = 
0, we have e^+^ = e^e^g-I^-^l/^, thus Eq.® can be 
put into the following form 

5-1 (r) |0(;3)) = sechfeech^/Vexpl^^^^atat 



tanh r 



\ cosh? 

(5t2tanh2 0-at2)| |00). (9) 



Next, we shall use Eq.® to derive the average of 
operator e-*" " under the suqeezed thermal vacuum 
state 5*1 (r) |0(/3)), which is a bridge for our whole 
calculations. Notice e^/'^^^^a'^e-f^'^'''''' = ah^^^ and 

g-//2a a^g//2a a _ fjg//2 gQ .^^g J^avC 



3//2ata^^ (r) |0(;3)) = sech6'sechi/ V 



exp 



f tanh 9 



\ cosh 



t;;t,//2 



-a' a' e 



tanhr 



tanh^ 9 



|06), 
(10) 



which leads to 



(0(/3)|5|(r) e^^'-^Si (r) |0(/3)) 

■ /2 tanh 9 



2 / ~ I I * 

— sech ^sechr (00| exp < - 



coshr 



tanh? 



X exp 



f tanh 9 



\ cosh 
tanhr 



h^^tanh^9-a^^e^) \ |00) 



= [Ce'^f - 2Bef + A] 
where we have set 



^1/2 



(11) 



III. SUQEEZED THERMAL VACUUM STATE 

To realize our purpose, we first introduce the squeezed 
thermal vacuum state, defined as (r) |0(/3)) , where 
Si (r) = exp[r(a^ — a^^)/2] is the single-mode squeez- 
ing operator for the real mode with r being squeezing 
parameter. Note that Eq.Q and the Baker-Hausdorf 
lemma 

(7) 
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C = cosh^ 



cosh 9 -I- cosh 29 sinh r 
^ ' (2??c + 1) cosh^ r, 

cosh 29 cosh^ r 
(2r?c + 1) sinh^ r. 
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sinh^ 9 cosh^ 9 
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5*1 (r) Sl (r) — coshr -t- asinhr, 
then we get 

Si (r) |0(/3)) = 5i (r) sech6'exp [a^'d'' tanh6l] |o6) 

= sech6'sech^^^r exp [(a^ coshr -I- asinhr) tanh 



and we have used the completness relation of coherent 
state / (Pzcfz\zz) (zz| /tt^ = 1, here \z) and \z) is the 
coherent state in real and fictitious modes, respectively, 
and the formula 1251 



J ^ exp (C \zf +^z + ?;z* + fz^ + gz*^^ 



X exp 



at2 ■ 

tanh r 

2 



1 00) 



(8) 
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whose convergent condition is Re(C ± / ± g) < 0, Re(C^ — 
4/5)/(C± /iff) < 0. Eq.([ni) is very important for later 
calculation of photon-number distribution (PND) and 
normalization of photon-added (-subtracted) squeezed 
thermal states (PASTS, PSSTS). 

It is interesting to notice that the standard generating 
function of Legendre polynomials [26| is given by 



Note the operator identity e'^"^" = e"'^:cxp[(l - 
e^'^)a^a]:, we see 

^-a"V"-;e-^°;=^-^ , (18) 

m=o ' ^ ^ 



1 



Vl - 2xt + t"^ 



(14) 



m— 



thus comparing Eq. lfTTI) with Eq. (fT4l) we find 



Ja^a 



^e^f _ 2^ef + 1 
A A 



-1/2 



^-1/2 



m=0 



(15) 



which indicates that the average value of e-*" " under 
squeezed thermal state (STS) is just the generating func- 
tion of Legendre polynomials, a remarkable result. Next, 
we shall examine the normalizations and PNDs of PASTS 
and PSSTS by using Eqs.dH]) and 



IV. NORMALIZATION AND PND OF PASTS 

The m-photon-added scheme, denoted by the mapping 
p — > a^™pa''", was first proposed by Agarwal and Tara 2] . 
Here, we introduce the PASTS. Theoretically, the PASTS 
can be obtained by repeatedly operating the photon cre- 
ation operator on a STS, so its density operator is 
given by 



Pad = C-la^"'SipthSla"', 



(16) 



where m is the added photon number (a non-negative 
integer), C~,^„ is the normalization constant to be deter- 
mined. 



A. Normalization of PASTS 

To fully describe a quantum state, its normalization is 
usually necessary. Next, we shall employ the fact ([5]) and 
([TT|). ([T5|) to realize our aim. According to the normal- 
ization condition txpad = 1 and the TFD, we have 

Ca,m = tr [at™^,p,^5|a"" 

= Tr [at™5i|0(/3))(0(/3)|5iV 

= (0(/3)|5iW™5i|0(/3)), 



(17) 



which implies that the calculation of normaliation factor 
C'a.m is converted to a matrix element after introducing 
the thermal vacuum state |0(/3)). 



where the symbol : : denotes antinormally ordering. Thus 
using Eqs.dn]), ([HI) and {e^ ^ j^) we have 



E ^^^^^rn = {m\ sU'^'-'-'^s, \m) 



m=0 



where we have set 



[At^ - 2Dt + 1] 



-1/2 



(19) 



D — cosh^ 9 cosh 2r ~ sinh^ r 
= ric cosh 2r + cosh^ r. 



(20) 



Comparing Eg. dTOl) with Eg. dMl) . and taking t' 
At, we obtain 



m=0 



\t''^ -2DIVAt' + \ 



-1/2 



= E ^™ (^/^) (21) 

thus the normalization constant of PASTS is given by 



(22) 



which is identical with the result in Ref.[28|. It is 
noted that, for the case of no-photon- addition with 
TO = 0, Cq.o = 1 as expected. Under the case of 
TO-photon-addition thermal state (no squeezing) with 
D ^ Uc + 1, A ^ {uc + lf , and P,„ (1) = 1, then 
Ca,m — "T'! {n-c + 1)™ ■ The same result as Eq.(32) found 
in Ref . [29| . In addition, when r = corresponding 
to photon-added thermal state, Eg. ip^ just reduces to 

Ca,m = TO!cOSh2™6' 



B. PND of PASTS 

The photon-number distribution (PND) is a key char- 
acteristic of every optical field. The PND, i.e., the prob- 
ability of finding n photons in a quantum state described 
by the density operator p, is V{n) = tr [|n) (n| p]. In 
a similar spirit of deriving Eg. (1221) . noting a™ |n) = 
v/n!/(n-m)!|n-m) and |n) = a'f''/V^.\0) , the PND 
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of the PASTS can be calculated as 

^a(n) = C-,i„tr \\n) {n\ a^'^SipthSla" 



C„ i,Tr 



\n) (n|at"^i|0(/3))(0(/3)|5ja" 



^^(0(/3)|5j|0 {l\S,\m) 
n\C-l . 



(0(/3)|5iV|0)(0|a'^i|0(/3)) (23) 



which leads to 



1=0 



oo 1 

— r' 



a'a e : 



5i |0(/3)) 



m)\si 



t a 111 



"^1 |o(/?)) , 



(24) 



where I = n — m and the vacuum projector op- 
" : and the operator identity 
a] : are used. 

Using Ea. ([TT|) again {e^ — > t) and comparing Eq. ([M)) 
with Eq. (|14p we see 



erator |0) (0| e"*^ ' 
gAata ^. exp[(e^ - l)at 



II 



!=0 



-rT - 2— T + 1 



-1/2 



which leads to the PND of PASTS 

n\C-UCIAt--^l^ 



Vain) 



[n - m)!vA 



'AC 



(25) 



(26) 



a Legendre polynomial with a condition n ^ m which 
implies that the photon-number (n) involved in PASTS 
is always no-less than the photon-number (m) operated 
on the STS, and there is no photon distribution when 
n < m). It is obvious that when m — corresponding 
to the STS, then the PND of STS is also a Legendre 
distribution Isoll. 



V. NORMALIZATION AND PND OF PSSTS 

Next, we turn our attention to discussing the photon- 
subtracted squeezed thermal state (PSSTS), defined as 



Psb = C^l^a"'SipthS[a 



t„tm 



(27) 



where m is the subtracted photon number (a non- 
negative integer), and Cg.m is a normalized constant. 



In a similar way to deriving Eq. (j22p . we have 

C,,™ = (0(/3) I Sy^^a^Si \0m, (28) 

so employing e'^°*° —: exp[(e'^ — l)a^a]: and Eq. lfTTj) 
(e^ ^ 1 + r) we see 



^ m 



m=0 



= [Ct^ - 2Et + 1] 



-1/2 



(29) 



where 



E = cosh 2r cosh^ 9 — cosh^ r 



- [(2nc + l)cosh2r - 1] 



Comparing Eq.([2H]) with Eq.([Tl|) yields 
= m!C™/2P™ (e/VC 



(30) 



(31) 

which is the normalization factor of PSSTS. When 
r — corresponding to photon-subtracted thermal state, 
Eq.dHl]) just reduces to C^,™ = ml sinh^" 9 

Using the same procession as obtaining Ea. (PB|) . the 
PND of PSSTS is given by 

-Psin) = C-i (0(/3)| 5lat" \n) {n\ a^S, |0(/3)) 

= ^C^l (0(/3)| S\ : at'"+"a™+"e-'^''^ : 5i |0(/3)) , 

(32) 

SO [k — m -\- n) 

oo ^ 

^n!C,,„P,(7^) = (0(/3)| S^^^'^-S, |0(/3)) 



fe=0 



= R.H.S of (EH), 
whcih leads to the PND of PSSTS 



(33) 



(34) 

a Legendre polynomial, which is same as the result of 
Ref.fsoJ. 



VI. CONCLUSION 

In this paper, we present a new concise approach for 
normalizing m-photon-added (-subtracted) STS and for 
derivingthe PNDs, which improve the method used in 
Refs. [2^, [S^j. That is to say, using the thermal field 
dynamics theory, we convert the thermal state to a pure 
state in doubled Fock space in which the calculations 
of ensemble averages under a mixed state p, {A) = 
tr {Ap) /tr [p) is replaced by an equivalent expectation 
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values with a pure state |0(/3)), i.e., {A) = (0(/3)| A |0(/3)). 

It is shown that the average value of e-'" " under STS is 
just the generating function of Legendre polynomials, a 
remarkable result. Based on this point, the normaliza- 
tion and PNDs of m-photon added (or subtracted) STS 
are easily obtained as the Legendre polynomials. The 

generating function of the Legendre polynomials and the 

f t 

average value of e-''^ under STS are used in the whole 
calculation. 
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